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MULTI- ACCESS  STRATEGIES  FOR  VOICE/DATA  INTEGRATION 
IN  HETEROGENEOUS  MIXED-MEDIA  PACKET  RADIO  NETWORKS 

1.  INTRODUCTION 

In  a  variety  of  commercial  and  military  communication  networks,  it  is  possible  that 
the  network  links  employ  different  communication  media  with  markedly  different  char¬ 
acteristics.  The  different  characteristics  could  be  different  frequency  bands,  data  rates 
(channel  capacities),  error  rates,  delays,  and  network  connectivity.  We  call  these  networks 
heterogeneous  mixed-media  networks,  while  we  use  the  terminology  of  multi-media 
networks  to  refer  to  networks  supporting  different  types  of  traffic,  such  as  voice  and  data. 
The  networks  studied  in  this  report  actually  fall  into  both  of  these  classes. 

As  examples  of  such  heterogeneous  mixed-media  networks  consider  radio  net¬ 
works,  which  use  the  HF  (both  groundwave  and  skywave),  UHF,  and  EHF  SATCOM 
frequency  bands.  In  these  networks,  some  of  the  links  are  terrestrial  links  (not  necessarily 
line-of-sight),  while  others  are  satellite  links.  The  scenario  for  the  ground  component  of 
such  networks  may  be  of  terrestrial  based  packet  radio  integrated  networks  that  support 
mobile  terminals.  They  can  also  be  a  subnet  of  a  larger  system  of  interconnected  integrated 
networks  that  provide  radio  interface  for  some  users  of  a  backbone  packet-switched  net¬ 
work.  The  satellite  component  may  involve  a  single  satellite  over  the  geographical  area  of 
the  terrestrial  subnetwork  or  even  a  network  of  multiple  satellites  covering  a  much  br  ader 
geographical  area  with  many  distinct  terrestrial  networks. 

Mixed-media  networks  are  important,  because,  if  suitable  multiple-access  and  routing 
algorithms  are  devised,  the  overall  network  efficiency  can  be  considerably  increased  by  using 
all  available  resources  (media,  e.g.,  communication  links  over  different  frequency  bands). 
Furthermore,  for  military  applications  and  under  stressed  conditions,  the  combined  use  of 
all  media  may  provide  connectivity,  whereas  use  of  a  single  medium  may  not. 

In  the  paper  by  D.  Huynh,  H.  Kobayashi,  and  F.  F.  Kuo  ([1]),  a  mixed-media  network 
was  considered  consisting  of  a  ground  subnet  and  a  satellite  subnet.  The  ground  subnet 
was  connected  with  wired  lines  and  used  standard  store-and-forward  packet  switching.  The 
satellite  subnet  used  contention-based  channel  access  (ALOHA).  The  retransmissions  of 
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the  packets  that  collided  in  the  satellite  channel  were  split  between  the  ground  links  and 
the  satellite  link.  Only  data  packets  were  transmitted  through  the  network  of  [1]. 

The  other  issue  of  interest  to  this  report  is  multi-media  integration  over  a  common 
channel.  This  is  motivated  by  evolution  of  communication  networks  toward  integrated 
services  digital  networks  (ISDN)  to  accommodate  random  demands  for  service  from  a 
population  of  heterogeneous  users  with  a  variety  of  data  types  (such  as  interactive  data, 
digital  voice,  video).  These  traffic  types  (in  this  report  we  restrict  attention  to  voice 
and  data)  have  different  characteristics;  voice  can  tolerate  higher  error  rates  but  requires 
constant  delays  with  no  variability,  whereas  data  can  tolerate  longer  delays  but  require 
lower  error  rates. 

Integration  of  different  traffic  types  over  a  common  channel  requires  a  method  to 
determine  how  the  users  should  coordinate  their  transmissions  to  avoid  the  destructive 
interference  which  occurs  when  the  traffic  load  over  the  channel  increases.  In  other  words, 
a  channel-access  scheme  has  to  be  followed  by  the  terminals  of  possibly  different  traffic 
types  in  order  to  make  efficient  use  of  the  channel.  In  a  multi-media  environment  hybrids 
of  circuit  switching  and  packet  switching  techniques  have  been  recommended  (see  [2]- [4] ) 
which  accommodate  both  the  need  of  voice  traffic  for  real-time  delivery  and  the  required 
low  packet  error  probabilities  of  data  traffic. 

In  the  recent  work  of  [5]  fixed-  and  movable-boundary  channel-access  schemes  for 
voice/data  integration  in  wireless  networks  were  introduced  and  their  performance  ana¬ 
lyzed.  Of  particular  interest  were  the  movable- boundary  schemes  of  [5].  which  allocate  in  a 
dynamic  manner  time  slots  to  voice  or  data  traffic  according  to  the  overall  network  traffic 
conditions.  Moreover,  in  our  recent  work  of  [Gj  multiple-access  schemes  for  voice/data  inte¬ 
gration  in  spread-spectrum  packet  radio  networks  were  introduced  and  their  performance 
evaluated.  The  introduction  of  movable-boundary  schemes  in  the  code  domain  was  the 
gist  of  [G].  The  work  of  both  [5]  and  [G]  pertained  to  networks  with  a  single  medium. 

In  this  report,  we  consider  a  packet  radio  network  with  ground  subnet  using  radio 
links  and  a  satellite  subnet.  For  example,  HF  links  may  be  used  for  the  ground  subnet  and 
EHF  SATCOM  links  for  the  satellite  subnet.  Both  voice  and  data  packets  are  transmitted 
through  the  network.  The  terrestrial  (or  ground)  channel  is  characterized  by  a  bandwidth 
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W\,  delay  rj,  and  packet  error  rate  Pe,i-  The  corresponding  features  of  the  satellite 
(SATCOM)  channel  are  W2,  r2,  and  Pe,2,  respectively.  It  is  assumed  that  rj  <  r2  and 
Pe,  1  >  Pe,2 •  These  comparisons  are  typical  of  the  relationships  of  the  two  channels  in 
practical  situations. 

In  this  report  efficient  multiple-access  schemes  are  devised  which  take  advantage  of 
the  characteristics  of  the  terrestrial  and  satellite  channels  above  and  accommodate  the 
requirements  of  the  different  traffic  types  (voice  and  data).  In  accomplishing  this,  our 
report  extends  the  work  of  (1]  for  purely  data  mixed-media  networks  to  radio  mixed- 
media  networks  with  two  traffic  types  (voice  and  data),  our  work  of  [6]  for  single-medium 
networks  to  mixed-media  networks,  and,  to  a  certain  extent  the  work  of  [5]  for  single¬ 
medium  networks  to  the  mixed-media  case. 

In  the  terrestrial  radio  network  Code-Division  Multiple- Access  (CDMA)  is  used.  This 
is  motivated  by  the  fact  that  CDMA  can  support  a  good  number  of  asynchronous  (at 
the  frequency- hop  or  dwell-time  level)  simultaneous  users,  it  mitigates  the  effects  of  fad¬ 
ing  and  narrowband  interference,  it  combats  jamming,  and  it  provides  selective  address- 
ing/reception  capability.  Moreover,  the  performance  of  CDMA  schemes  degrades  grace¬ 
fully  as  the  interference  level  (from  any  source)  increases;  this  can  be  used  efficiently  for 
voice  traffic  given  the  graceful  degradation  of  voice  quality  with  increasing  in-band  (or 
co-channel)  interference.  Similar  feedback  and  voice/data  integration  schemes  as  those 
developed  in  [6]  (termed  schemes  A,  B.  and  C  in  that  report)  can  be  used  for  the  ground 
channel. 

On  the  other  hand,  a  framed  protocol  employing  reservations  for  the  voice  traffic  and 
ALOHA  with  retransmission  control  for  the  data  traffic  is  used  by  the  satellite  channel: 
this  is  a  version  of  the  framed  ALOHA  protocol  (see  [7]- [8]).  modified  to  accommodate 
both  voice  and  data  traffic. 

The  choice  of  channels  and  basic  transmission  policies  for  the  two  traffic  types  is 
described  and  justified  in  the  following.  The  voice  traffic  uses  primarily  the  terrestrial 
subnetwork  and.  if  necessary,  the  satellite  subnetwork  as  well.  This  is  motivated  by  the 
requirement  for  short  delay  times  and  low  variability  of  delay  of  the  voice  calls.  Also 
a  higher  level  of  interference  is  present  in  the  terrestrial  channel  making  it  less  reliable. 
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However,  since  voice  can  tolerate  higher  levels  of  interference  than  data,  the  selection  of 
the  terrestrial  network  as  the  prime  carrier  of  voice  traffic  is  well  justified. 

Due  to  the  more  stringent  requirements  for  smaller  packet  error  probability  and  the 
tolerance  of  longer  delays  of  the  data  users,  data  traffic  is  first  transmitted  through  the  more 
reliable  satellite  channel  and,  in  the  event  of  packet  collisions,  some  of  the  retransmitted 
data  packets  are  sent  through  the  terrestrial  subnetwork,  while  the  remaining  packets  are 
sent  through  the  satellite  channel.  This  splitting  is  performed  dynamically  based  on  the 
current  state  of  voice  and  data  traffic  in  the  two  channels. 

The  mixed-media  network  of  interest  to  this  report  is  illustrated  in  Figure  1  with 
examples  of  different  types  of  links. 

In  our  report  we  study  the  tradeoffs  among  the  following  performance  measures: 
throughput  of  data  and  voice  packets,  delay  of  data  packets,  and  voice  call  blocking  prob¬ 
ability.  The  analysis  is  in  principle  applicable  to  several  different  strategies  for  upgrading 
the  retransmission  probabilities  and  the  splitting  probabilities:  however,  for  the  sake  of 
clarity  in  the  presentation  and  simplicity  in  generation  of  the  numerical  results  only  one 
such  strategy  is  analyzed  and  and  evaluated. 

This  report  is  organized  as  follows.  Section  2  presents  the  system  model  and  the 
detailed  protocol  operation  for  the  ground  and  satellite  subnetworks.  In  Section  3.  we 
develop  a  Markovian  model  for  voice  traffic  analysis.  Section  4  follows  with  the  model  and 
the  analysis  for  data  traffic.  In  Section  5  numerical  results  about  the  system  performance 
are  presented. 

2.  SYSTEM  MODEL 

We  consider  a  packet  radio  network  with  a  population  of  Md  data  users  and  M, 
voice  users.  The  network  consist  of  a  ground  subnet  (G-subnet ).  using  ground  radio  links, 
and  a  satellite  subnet  (S-subnet).  The  time  axis  of  the  satellite  channel  is  grouped  into 
frames  of  length  L  slots  each.  The  frame  length  should  be  at  least  as  large  as  the  round 
trip  propagation  delay  of  the  satellite  link;  this  is  elaborated  upon  in  Section  2.1  below. 
Time  over  the  ground  channel  is  also  measured  in  frames  but  there  is  only  one  time  slot 
in  each  frame.  The  basic  time  unit  is  the  slot  which  represents  the  time  allowed  for  the 
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transmission  of  a  single  data  or  voice  packet.  The  frames  of  the  ground  channel  and  the 
frames  of  the  satellite  channel  are  perfectly  aligned.  It  is  assumed  that  the  bandwidth 
of  the  satellite  channel  is  equal  to  L  times  the  bandwidth  of  the  ground  channel  for  the 
information  signals.  Therefore,  the  information  rate  over  the  satellite  channel  is  L  times 
the  information  rate  over  the  ground  channel  and,  since  each  user  transmitting  over  the 
satellite  channel  uses  only  one  slot  in  a  frame,  the  same  amount  of  information  per  user  is 
transmitted  during  one  frame  of  the  satellite  channel  and  one  frame  of  the  ground  channel. 
Figures  2  and  3  illustrate  the  protocol  operation  over  the  satellite  and  ground  subnets, 
respectively,  over  three  consecutive  frames.  The  detailed  operation  is  described  below. 

2.1.  Channel  Access  for  Voice  Traffic 

Code-division  multiple-accessing  (CDMA)  is  used  over  the  ground  subnet  to  accom¬ 
modate  several  voice  calls  simultaneously.  In  our  model,  packet  transmissions  start  at 
common  clock  instances  and  packets  have  constant  length  of  N  symbols.  Each  user  (data 
or  voice  node)  employs  a  random  frequency  hopping  pattern  for  the  transmission  of  its 
packets.  The  frequency  spectrum  is  divided  into  q  frequency  bins  and  each  symbol  is 
transmitted  at  a  frequency  chosen  from  the  q  frequencies  with  equal  probability,  inde¬ 
pendently  of  the  frequencies  chosen  for  other  symbols.  We  assume  that  some  form  of 
Forward  Error  Control  (FEC)  coding  is  employed,  for  which  up  to  v  symbol  errors  can 
be  corrected.  We  also  assume  that  a  packet  consists  of  exactly  one  codeword  with  .V 
code  symbols.  Therefore,  a  packet  is  declared  successfully  transmitted,  if  at  most  u  sym¬ 
bol  errors  occur.  In  particular,  if  Reed-Solomon  (A'.  A’)  codes  with  erasures-correction 
decoding  are  employed,  i /  —  A  —  K  erasures  can  be  corrected.  Typically.  M-ary  FSK 
(Frequency-Shift-Keying)  data  modulation  with  M  =  A'  and  noncoherent  demodulation 
is  employed  in  frequency-hopped  spread-spectrum  (FH/SS)  systems.  I\  M- ary  symbols 
are  encoded  on  A’  Reed-Solomon  code  symbols  and  each  code  symbol  is  transmitted  on  a 
single  frequency  hop.  Each  voice  packet  of  a  typical  user  occupies  exactly  one  frame  over 
the  ground  channel. 

We  define  the  Multiple  Access  Capability  (MAC)  index  A"’  as  the  number  of  voice 
users  that  can  be  accommodated  simultaneously  such  that  the  expected  packet  error  prob- 
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ability  of  voice  traffic  remains  below  a  specified  threshold.  Specifically,  the  MAC  indices 
K%  and  K9d  (for  the  data)  are  defined  by 


PE(k)  <  P'e 

V  k  <  K% 

(2.1) 

PE(k)  <  PdE 

V  k  <  K9 

(2.2) 

where  Pg  and  Pg  are  the  maximum  tolerable  voice  and  data  packet  error  probabilities, 
respectively,  and  Pe( k)  is  the  packet  error  probability  in  the  presence  of  k  simultaneous 
packet  transmissions,  where  k  includes  both  voice  and  data  users.  In  practice,  P£  >  P£  and 
therefore  K9  >  Kjj.  If  the  total  number  of  simultaneous  users  is  k  <  A'J,  then  all  k  data  or 
voice  packets  are  received  with  acceptable  error  probability;  if  I\d  <  k  <  A'*?,  then  among 
the  k  packets  the  voice  packets  are  received  with  acceptable  error  probability,  whereas 
the  data  packets  have  higher  than  the  acceptable  error  probability;  finally,  if  A®  <  k. 
then  all  voice  or  data  packets  are  received  with  unacceptable  error  probabilities.  The 
model  described  above  is  referred  to  as  the  threshold  model  for  other-user  interference 
in  CDMA  systems. 

In  our  analysis  of  data  traffic  over  the  ground  subnet  we  also  use  the  graceful  degra¬ 
dation  model  for  CDMA  multiple-reception  of  data  packets.  According  to  this  model, 
there  is  a  non-zero  probability  of  correct  reception  for  any  arbitrary  number  of  pack¬ 
ets  which  depends  on  the  total  number  of  simultaneously  transmitted  packets  (data  and 
voice)  over  the  G-subnet.  This  model  is  described  in  detail  in  Section  4.2. 

Notice  that  this  discussion  assumes  the  same  code  rate  for  both  voice  and  data  packets. 
It  could  be  extended  to  systems  with  different  code  rates  for  the  data  and  voice  packets. 

Up  to  A  *  number  of  slots  per  frame  may  be  allocated  to  the  voice  traffic  on  a  reser¬ 
vation  basis  over  the  satellite  subnet.  One  voice  packet  of  a  typical  user  occupies  exactly 
one  slot  per  frame  over  the  satellite  channel.  The  first  slot  of  each  frame  over  the  satellite 
subnet  is  used  by  the  voice  nodes  to  reserve  a  slot  for  packet  transmission  in  the  next 
frame.  We  assume  that  the  transponder  employed  for  the  satellite  subnet  is  located  on  a 
geosynchronous  satellite  which  results  in  a  round-trip  propagation  delay  of  0.27  seconds. 
Therefore,  the  frame  length  must  be  as  large  as  the  round-trip  propagation  delay,  so  that 
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the  reservation  information  will  be  available  to  the  users  before  the  start  of  the  next  frame. 
The  round-trip  propagation  delay  is  of  the  order  of  10-13  slots  for  typical  satellite  networks. 

Each  idle  voice  node  becomes  active  in  a  frame  with  probability  P„  and  attempts 
transmission  at  the  beginning  of  the  next  frame.  Each  voice  call  generates  a  random 
number  of  packets  geometrically  distributed  with  parameter  p.  This  means,  that  the 
probability  that  a  voice  call  lasts  for  /  frames  (transmitting  one  packet  per  slot  per  frame 
over  the  satellite  channel  or  one  packet  per  frame  over  the  ground  channel)  is 

P{ voice  call  duration  =  /  packets}  =  (1  —  p)l~lp  for  Z  =  1,2,---  (2.3) 

The  value  l  =  0  was  excluded  from  (2.3),  because  it  is  assumed  that  a  call  that  just  became 
active  can  not  terminate  before  it  lasts  for  at  least  one  packet  (frame).  Thus,  the  average 
duration  of  a  voice  call  is 

V/(l  -p),-1p=  -  >  1  (2.4) 

fcl  P 

which  is  larger  than  one  packet. 

It  is  assumed  that  some  information  about  the  number  of  established  voice  calls  in 
progress  is  available  to  the  nodes  of  the  network  at  the  beginning  of  each  frame.  This 
information  may  be  obtained  by  observing  (sensing)  the  channel  activities  during  each 
frame.  Several  sensing  schemes  are  possible  that  provide  different  amounts  of  information 
about  the  state  of  the  channel. 

Among  the  many  possibilities,  three  different  feedback  and  control  schemes  were  con¬ 
sidered  in  [5]  as  the  most  representative  and  meaningful.  Scheme  A  characterized  the 
situation  where  the  nodes  of  the  network  observe  the  channel  activities  in  each  slot  and 
obtain  some  feedback  information  based  on  the  number  of  combined  active  voice  and  data 
terminals.  In  Scheme  B.  it  was  assumed  that  a  central  controller  provides  the  nodes  of 
the  network  with  the  exact  number  of  established  voice  calls  and  data  nodes  in  backlog. 
Scheme  C  is  a  variation  of  Scheme  B.  where  the  voice  terminals  do  not  use  the  information 
about  the  data  backlog  to  regulate  their  access  to  the  channel.  Variations  of  Schemes  A. 
B.  and  C  were  also  considered  according  to  which  data  and  voice  share  the  multiple-access 
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capability  of  the  CDMA  channel  in  a  movable-boundary  fashion  with  varying  degrees  of 
fairness  in  the  initial  and  time-varying  allocation  of  the  MAC  between  them. 

In  this  report  we  assume  that  channel  sensing  (feedback)  makes  available  (i)  the 
number  of  ongoing  voice  calls  over  both  the  G-subnet  and  the  S-subnet  to  all  voice  users, 
and  (ii)  the  data  backlog  over  the  G-subnet,  the  number  of  colliding  data  packets  over 
the  S-subnet,  and  the  number  of  ongoing  voice  calls  over  the  S-subnet  and  the  G-subnet 
to  all  data  users.  The  voice  users  do  not  use  information  about  the  data  traffic.  This  is 
a  variation  of  Scheme  C  of  [5]  modified  to  be  suitable  in  mixed-media  (ground/satellite 
networks)  environment. 

Voice-call  blocking  is  used  as  a  mechanism  to  control  the  flow  of  voice  traffic.  At  the 
beginning  of  each  frame,  the  voice  nodes  attempting  to  establish  new  calls  observe  the  state 
information  of  the  ground  channel.  If  the  number  of  established  voice  calls  over  the  ground 
subnet  is  equal  to  the  MAC  index  A'*?,  the  voice  nodes  observe  the  state  information  of  the 
satellite  channel.  If  the  number  of  established  calls  over  the  satellite  channel  equals  the 
maximum  capacity  allocated  to  the  voice  traffic  over  the  satellite  subnet  (i.e..  A'®),  then 
the  new  calls  will  be  blocked  and  cleared.  Otherwise,  up  to  A'*  number  of  calls  may  be 
established.  With  this  access  protocol  for  voice  traffic,  once  a  voice  call  is  established,  the 
voice  packets  either  experience  no  delays  (i.e.,  voice  calls  over  the  ground  subnet),  or  have 
a  constant  delay  of  L  slots  (i.e.,  voice  calls  over  the  satellite  subnet). 

2.2.  Channel  Access  for  Data  Traffic 

Due  to  the  high  reliability  of  the  satellite  channel,  the  data  traffic  primarily  goes 
through  the  satellite  subnet.  Actually,  all  data  nodes  initially  transmit  over  the  satel¬ 
lite  channel.  At  the  beginning  of  each  frame,  the  data  nodes  have  information  about  the 
number  of  established  voice  calls  in  that  frame.  This  information  determines  the  number 
of  slots  available  for  data  ■  acket  transmission  in  that  frame.  We  assume  that  each  data 
terminal  generates  a  new  packet  (independently  of  previous  packets  and  other  terminals) 
with  probability  P*  in  each  frame  and  remains  idle  with  probability  1  —  Pd-  It  is  thus 
possible  that  a  data  node  generates  a  new  packet  in  a  frame  although  some  of  its  previ¬ 
ously  generated  packets  have  not  yet  been  transmitted  successfully.  It  is  assumed  that  all 
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data  nodes  have  at  their  disposal  infinite  size  buffers  (moderately  large  size  buffers  will 
suffice  in  most  practical  situations)  where  they  store  the  backlogged  packets  waiting  for 
retransmission  and  the  newly  arrived  packets  waiting  for  their  first  transmission. 

The  framed  ALOHA  protocol  is  used  by  the  data  nodes  to  access  the  satellite  channel 
[7]- [8].  All  the  new  data  packet  arrivals  during  a  frame  will  be  transmitted  during  that 
frame.  The  result  of  transmissions  of  data  packets  will  be  known  by  the  end  of  the  next 
frame.  The  data  terminals  that  were  involved  in  a  collision  event  may  retransmit  the 
collided  packets  either  through  the  satellite  subnet  or  the  ground  subnet.  The  splitting 
of  the  retransmission  traffic  between  the  satellite  and  ground  subnet  may  be  controlled 
dynamically  at  the  beginning  of  each  frame.  The  data  packets  that  join  the  backlog 
of  data  packets  over  the  ground  subnet  are  transmitted  employing  an  ALOHA  scheme 
with  retransmission  control.  Since  CDMA  is  used  over  the  G-subnet  that  can  support  a 
maximum  of  Kd  simultaneous  data  and  voice  packet  transmissions  at  acceptable  packet 
error  probabilities,  the  retransmission  probabilities  for  each  frame  depend  on  I\9d  and  the 
number  of  ongoing  voice  calls  over  the  ground  network  during  this  frame. 

3.  VOICE  TRAFFIC  ANALYSIS 

We  begin  with  the  introduction  of  random  processes  that  model  the  voice  traffic. 

Let 

X,  g  —  Number  of  voice  calls  in  progress  in  the  G-subnet  at  the  beginning  of  frame  ; 

N't  *  =  Number  of  voice  calls  in  progress  ;n  the  S-subnet  at  the  beginning  of  frame  t 

.4)'  =  Number  of  new  voice  call  arrivals  in  frame  i 

E 9  =  Number  of  calls  accepted  by  the  G-subnet  in  frame  i  and  transmitted  in  frame 

i  +  1 

Xf  -  Number  of  calls  terminated  in  frame  i  over  the  G-subnet 

E'  =  Number  of  calls  accepted  by  the  S-subnet  in  frame  t  and  transmitted  in  frame 

i  +  1 

A’,"  =  Number  of  calls  terminated  in  frame  i  over  the  S-subnet 
Figure  4  depicts  all  components  of  voice  traffic  over  the  ground  and  satellite  subnets. 
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The  evolution  of  {X”9}  and  {N-13}  is  described  by 


=  N?9  +  E?  -  X? 

Nf+y  =  N™  +  E3  -  X*  (3.1) 


where 


'  0 

xr  >  k9 

Ef=< 

-4r_,  <  k*  -  xr,  xr  <  x> 

(3.2  a) 

i 

[  K?  -  X*9  Al’_x  >  A*®  -  A7S,  X’9  <  K? 

r  av 

■Hi-i 

.4j’_,  <  k3  -  xr ,  xr  >  k9 

K3  -  X r 

a*_i  >  i<t  -  Nr<  ~K9  >  k9 

e:=< 

0 

,4*’_,  <  x9  -  xr ,  xr  <  k9 

K-x-iKl- 

-K9) 

o  <  .4;_,  -  (K?  -  xr)  <  x*  -  xr ,  xr  < 

A'® 

.  k3v  -  xr 

au  -  (K9  -  xr)  >  k3  -  xr ,  Arr  <  a® 

(3.3a) 

or  equivalently, 

Ef  = 

r. 

.v;?  >  a'® 

(3.26) 

(  min(.4|. 

_i,  a'®  -  xr)  xr  <  *® 

niin(.4;_, .  A'*  -  A',"") 
min( i4J_,  -  Ef.A’J  -  A',M) 


X'  9  >  K9 

K9  <  A'f 


(3.3 b) 


Let  .V)' =(N-’a.  N*9).  From  the  above,  we  see  that  A’,1'*  and  X'9  are  dependent  on 
each  other,  but  the  joint  statistics  of  JVj'  can  be  shown  to  satisfy  the  Markovian  property. 
The  transition  probabilities  of  JV’’  denoted  by  P[X]'+i  =  (n,m)  j  A_J  =  ( ^' •> J ) ]  may 
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obtained  by  observing  that 


m?+1  =  (".">)  i  isr  =  (*.*)] 

=  P\Kil  =  Kt,  =  m  !  Nr  =  k,N‘a  =  ]} 

min{Mv—k  —  j,K’  —j,m  }  min{  Afv  —  k  —  j  —  /i  ,A‘*  —  k,n] 

53  E  p[  =  /i ,  p,j  =  h  i  Arr  =  = ;] 

/ !  =max  {0,m  —  j}  l;=mai(0,n-i) 

•  P[A7A  =  n,^+9,  =  m|£f  =  /i,£?  =  /2,  ^ 


The  first  term  in  the  sum  of  (3.4)  is  given  by 

P[JV,"  =  n.  A’,;5,  =  m  I  E?  =  f,  ,  £?  «  h  ,  AT,"  =  t.  A','*  =  j] 

=  PliV”;,  =  m  |  £f  =  A'.'*  =  jj  •  P[A\'V’,  =  n  |  £,■  =  J,.  A?*  =  t]  (3.5) 

Equation  (3.5)  holds  because  of  the  fact  that  when  conditioned  on  (Ef  —  1 1  ,  E*  =  1 2).  the 
number  of  accepted  calls  over  the  G-subnet  and  the  S-subnet,  and  on  ( A7”*  —  k,  A"**9  =  j ), 
the  number  of  ongoing  calls  at  the  beginning  of  frame  i ,  AT,'7i  and  A’t,+  1,  (the  number  of 
ongoing  calls  at  the  beginning  of  frame  i  +  1)  are  independent.  Recall  that  memory  less 
geometric  distributions  determine  the  duration  of  each  call. 

The  two  terms  in  the  right  hand  side  of  (3.5)  are  given  by 

P[A’2,  =  m  |  Ef  =  /,,  A7,"5  =  j]  =  b(j,m  -  /,,1  -  p)  (3.6) 


P[A7+S,  =  n  I  E;  =  /2,  ,V«  =  1']  =  b(k,n  -  /2. 1  -  p). 


where  the  binomial  distribution  is  defined  as 


MA'.n,p)=(‘Mp"(l-p)N  "  ;  0  <  n  <  A’. 


In  (3.6)  and  (3.7)  1  -  p  is  the  probability  that  an  active  call  in  the  S-subnet  during 
frame  1  continues  during  frame  i  +  1.  The  differences  m  —  li  and  n  —  /2  are  involved  in 
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(3.6)  and  (3.7),  respectively,  because  the  newly  accepted  calls  Ef  and  E 3  over  the  G-  and 
S-subnets  during  frame  i  are  actually  transmitted  in  frame  i  +  1  and  cannot  terminate 
before  transmitting  at  least  one  packet  according  to  the  model  of  (2.3). 

The  second  term  in  the  sum  of  (3.4)  can  be  obtained  from  equations  (3.2a)  and  (3.3a). 
As  seen  from  (3.2a)  and  (3.3a)  the  two  numbers  of  accepted  cedis  E9  and  E*  depend  not 
only  on  the  number  of  voice  calls  in  progress  A7**9  and  N?s  in  frame  i  but  through  the 
number  of  new  arrivals  they  depend  also  on  the  number  of  voice  calls  in  progress  Artv_?j 
and  N™]  in  frame  i  —  1.  To  capture  this  dependence  accurately  would  transform  the  Markov 
chain  of  (3.4)  from  a  first-order  to  a  second-order  one.  This  will  greatly  complicate  the 
analysis.  Instead,  we  make  the  following  well  justified  approximation.  Since  the  numbers 
of  voice  calls  in  progress  A7' 9  and  A7"3  change  very  slowly  from  frame  to  frame,  we  can 
substitute  A'9  and  A7,1’*  in  (3.2a)  and  (3.3a)  by  A’^,  and  A7’^,  respectively,  without 
substantial  loss  in  accuracy.  This  results  in 


P[  e9  -  h  .  e;  =  i2 1  av,j  =  k.  a7;9  =  j\  ^  p[  e9  =  /,  ,  e: 


&(A/„  -A’- J, /,,£,.) 
b(Mv  —  k  —  j,  1 1  -f  /  2 .  Pv) 

0 


1 1  <  A'9  -  j,  /2  =  o 
1 1  -  A'9  -  j,  /2  <  A'3  -  k 
h  =  A'9  —  j,h  =  A'3  -  k 

else  where 


(3.9) 


Once  the  transition  probability  matrix  of  AjJ’  is  obtained,  we  can  evaluate  the  station¬ 
ary  distribution  of  the  process  {II'’(A\  j )},  where  0  <  k  <  I\*  and  0  <  j  <  min  {A*9.  A/,.—  A-}: 
if  we  assume  that  A/,.  >  A'9  -f  A'3,  then  0  <  j  <  A’9.  The  voice  throughput  is  defined  as 
the  expected  number  of  voice  calls  in  progress  per  frame,  and  is  given  by 


r 7,.=  lim  £{A',''3  +  A*,"9}  =  V  +  j)Tl'(k.j) 

I  —  oc 


(3.10) 


<•  =  0  ;  =  0 


Moreover,  the  voice-call  blocking  probability  is  given  by 
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Pb 


E{ number  of  blocked  calls  per  frame} 
E{ number  of  arriving  calls  per  frame} 


lim  TT'Kt  +  KliT'Mv-n 
‘‘mi_OoLn=0  2^,1=  K‘  +  h  ‘  - 


lim,. 


-ooZ^  n=0 


n+1V-(K9  +  K!-n))P[A»_l=l\Nr  +  N»9  =  n}P[Nr  +  K9  =  n} 
E/=rn/pMr- 1  =  /|A^+A7s  =  n]P[A7*  +  A79  =  n] 


Enlo+A':E/=^n+^,-.+,[/-(^  +  A7  -n)]6(A/t.-n,/,P,)lim.,ocP[A7^A79^n 

Znlth\Mv  —  n)Pv  lim,— oo  P[N™  +  Kvl9  =  n] 


{3.11) 


The  above  equation  holds  under  the  ergodicity  assumption  for  the  Markov  Chain  ( A7*.  A79 ), 
so  that  the  time  averages  involved  in  the  practical  definition  of  blocking  probability  are 
equivalent  to  the  ensemble  averages  used  in  (3.11).  Also  we  used  the  well  justified  approx¬ 
imation 


P[.4;'_,  =  i\Nr  +  K9  =  n]  *  p[a;_j  =  /|A7J,  +  at_9,  =n)  =  hmv  -  nj,pv) 


according  to  the  discussion  preceding  (3.9).  A',?  -I-  A'*  represents  the  maximum  number 
of  voice  calls  that  can  be  accommodated  by  the  ground  and  satellite  subnets  (total  voice 
capacity  of  the  system).  According  to  (3.11),  blocking  occurs  when  the  number  of  new 
(arriving)  voice  calls  A}’_,  =  l  is  larger  than  the  system’s  voice  capacity  minus  the  number 
of  current  (ongoing)  voice  calls  A’,'9  +  A’,,s  =  n.  that  is.  when  l  >  A®  +  I\*  —  n ,  and  the 
numberof  blocked  calls  is  /-(A’-?  +  I\’—n).  The  stationary  distribution  of  P[A’7 9  +  A’’*  =  n] 
is  obtained  through  convolution  from  IT’^.j).  the  joint  distribution  of  ( ATIrs,  AT'9 ).  as 


min  {  A’f  ,n  } 

lim  P[A7'  +  A79  =  n]  =  V  IT(n  -  j.j  )  (3.12) 

1  —  OC  *' 

;  =  max  {0 ,n  —  i\*  } 


4.  DATA  TRAFFIC  ANALYSIS 

We  first  present  the  analysis  of  data  process  over  the  satellite  subnet. 

4.1.  Data  Traffic  over  the  5-Subnet 

At  the  beginning  of  frame  i.  all  data  nodes  have  information  about  the  number  of 
established  voice  calls  in  frame  i.  This  information  determines  the  number  of  slots  available 
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for  data  traffic  in  frame  i  given  by 


L,  =  L  -  1  -  A7* 


(4.1) 


Let 

l  =  Number  of  data  packet  transmissions  in  frame  i  over  S-subnet 
C,  =  Number  of  data  packets  involved  in  collision  in  frame  i  over  S-subnet 
S3  =  Number  of  successful  data  packet  transmissions  in  frame  i  over  S-subnet 
Nf3  =  Number  of  data  packet  retransmissions  in  frame  i  over  S-subnet 
=  Number  of  new  data  packet  arrivals  in  frame  i 
Figure  5  depicts  all  components  of  data  traffic  over  the  satellite  and  ground  subnets.  The 
necessary  notation  for  data  traffic  over  the  ground  subnet  is  provided  at  the  beginning  of 
Section  4.2. 

From  the  above  definitions,  we  see  that 


Y*  =  .4?  +  Xf3  =  S3  +  C , 


(4.2) 


It  is  assumed  that  all  the  data  nodes  that  were  involved  in  a  collision  event  in  frame  i  are 
able  to  estimate  the  number  of  packets  that  were  involved  in  collision  in  frame  i.  This 
information  becomes  available  at  the  end  of  frame  i  +  1  because  of  the  delay  (=  one  frame 
=  L  slots)  introduced  by  the  satellite.  Based  on  this  information,  a  decision  is  made  for 
each  collided  packet  at  the  beginning  of  frame  ?  +  2.  With  probability  ->,+2  it  is  transmitted 
over  the  S-subnet  sometime  during  frame  z  +  3  and  with  probability  1  —  7, +2  it  joins  the 
backlog  of  data  packets  over  the  G-subnet.  We  will  update  the  splitting  factor  •),.  based 
on  the  fraction  of  data  nodes  that  were  involved  in  a  collision  event.  This  rule  may  be 
given  by: 


(  1  C,<  aMd 


~U  +  2  = 


C’,  >  q  Md 


(4.3) 


where  a  >  0*  and  0  <  <  aMd ■  *  Clearly,  increasing  the  values  of  either  of  these 

parameters  results  in  an  increase  in  the  fraction  of  retransmissions  that  use  the  S-subnet. 

*  Note  that,  since  each  of  the  data  users  can  transmit  more  than  one  packet  per  frame, 
it  is  possible  to  have  C,  >  Md\  thus  a  >  1  may  be  appropriate  in  some  cases. 
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The  selection  of  the  splitting  probability  1  —  7, +2  above  is  not  optimal.  However, 
it  is  justified  by  the  fact  that,  as  the  number  of  colliding  data  packets  C,  increases,  it  is 
desirable  to  retransmit  a  larger  portion  of  these  packets  over  the  ground  subnet  so  that  the 
the  data  traffic  over  the  satellite  subnet  is  somewhat  relieved.  One  can  obtain  the  optimal 
values  of  the  parameters  o  and  /?  involved  in  (4.3)  through  a  numerical  search  for  the 
optimization  of  performance  criteria  such  as  the  data  throughput  and  average  data  delay. 

One  could  also  attempt  to  optimize  over  the  form  of  the  splitting  factor  7, +2,  but  this  is 
a  much  more  difficult  and  computationally  demanding  task  which  is  beyond  the  scope  of 
this  report. 

It  can  be  shown  that  the  evolution  of  the  data  process  over  the  S-subnet  can  be  de¬ 
scribed  by  the  two-dimensional  Markov  process  D*=(C,,  A7,1-* ).  The  transition  probability 
matrix  of  { D 3 }  may  be  obtained  by  observing  that 

P  Le;+2  =  (n,m)|£:  =  =  F  [C,+2  =  n|C,  =  k.N"  =  j.-V”,  =  m] 

'  p  [A’.+2  =  ”t\C,  =  =  >] 

=  P  [Ci+2  =  n|C,  .  k.  A’",  =  m]  •  P  [.V"2  =  m|.V,~  =  >j 

(4.4) 

This  decomposition  follows  from  the  fact  that  the  voice  process  N'3  is  independent  of 
the  data  process  and  the  number  of  data  packets  involved  in  a  collision  in  frame  t  4-  2  is 
independent  of  JVtvs  once  C,  and  Ar'+2  are  giyen. 

The  two-phase  transition  of  A\vs  can  be  obtained  from  the  2-step  transition  probability 
matrix  and  the  stationary  probability  distribution  of  Ay*'.  The  first  term  on  the  right  hand 
side  of  (4.4)  is  given  by 

4-  k 

P  [CI+2  =  n| C,  =  k\  A7+2  =  m]  =  P  [C,+2  =  «|C.  -  k.  av;%  =  m.r-;2  =  /] 

/=  n 

•  e  [v,;2  =  lie,  =  k.xz,  =  m] 

Mi  +  k 

=  Y.  p  =  "iA’.+2  =  ™.>r+2  =  1]  p [ir+s  =  'ic.  =  »■] 

/=n 
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for  n  <  Md  +  k ;  the  above  quantity  is  0  for  n  >  Md  +  k. 
The  second  term  in  the  sum  of  (4.5)  is  given  by 


(4.5) 


P  [Y‘+2  =  l\Ci  =  k]  =  Y,  P[K-2  =  a\C'  =  k]-P[N^2  =  l-a\C,  =  k]  (4.6) 

a  =  max(O.I- k) 

The  conditional  distributions  of  new  arrivals  and  the  retransmitted  packets  over  the 
S-subnet  are  given  by: 


P  [-4 ?+2  =  a\C,  =  k]  =  b(Md,a,Pd)  (4.7) 

P  [:V^2  =  /  -  a\Cx  =  k]  =  b(kJ-a.~,,+2)  (4.8) 

where  7,+2  is  given  by  (4.3).  Equation  (4.7)  is  justified  by  the  mechanism  of  data  packet 
generation  described  in  Section  2.2.  Since  the  Ct  data  terminals  with  colliding  packets 
can  generate  new  packets  in  frame  i  +  2  (at  the  time  they  find  out  about  the  status  of 
their  transmitted  packets  in  frame  t)  the  number  of  data  terminals  that  can  generate  new 
data  packets  is  Md  independently  of  the  value  of  C,.  On  the  other  hand  (4.S)  is  a  direct 
consequence  of  the  splitting  scheme  of  (4.3). 

The  first  term  in  the  sum  of  (4.5)  is  equivalent  to 

P  [cI+2  =  »|av+2  =  m,  r;+2  =  /]  =  p  [s;+2  =  /  -  n\N%2  =  m.  r,s+2  =  /] 

=  P  [S*+2  =  /  -  n|I”2  =  I  -  1  -  m.y;+2  =  /](4.9) 

where  the  conditional  distribution  of  S’  is  given  by: 


P  [S';  =  k\y,3  =  m.L,  -  j] 


,  i  l  l\  mln  \m'J  I  ,  .  \l ,  ■  ; \m  —  l 

(-1  rj.m\\  ^  (-1)0-0 _ 

jmk'  )  f-r  (/-  k)\{j  -mm  -  /)! 

J  '  l=k 


(4.10) 


Eq.  (4.10)  (taken  from  [9])  represents  the  probability  that  k  cells  are  occupied  with  one 
ball  (i.e,  the  successful  data  packet  transmissions  S’)  when  m  balls  (i.e..  the  transmitted 
data  packets  }’*)  are  thrown  into  j  cells  (i.e.,  the  number  of  available  slots  L,) 


1G 


An  important  comment  must  be  made  here  regarding  (4.10)  and  the  retransmission 
mechanism  of  the  framed  ALOHA  protocol  as  applied  to  our  mixed-media  network  model. 
Recall  that,  according  to  the  model  for  the  generation  of  new  data  packets  described  in 
Section  2.2,  it  is  possible  that  data  nodes  generate  new  packets  before  several  of  their  pre¬ 
viously  generated  packets  have  been  successfully  transmitted.  Thus,  several  data  packets 
belonging  to  the  same  data  node(s)  may  be  part  of  Nfa,  Ct,  S*,  and  V'*  of  (4.2).  Con¬ 
sequently,  a  number  of  data  packets  of  the  same  data  node  may  be  competing  for  the  L, 
[given  by  (4.1)]  available  slots  of  the  i-th  frame.  Under  these  conditions,  any  practical 
implementation  of  the  framed  ALOHA  protocol  will  use  some  mechanism  to  avoid  these 
suicidal  collisions  between  packets  belonging  to  the  same  data  node  by  spreading  the  trans¬ 
missions  of  these  packets  over  several  consecutive  frames  or  by  simply  choosing  different 
slots  within  (possibly)  the  same  frame  for  its  packets.  This  is  less  of  a  problem  for  larger 
values  of  Mj.  However,  the  analysis  of  such  a  scheme  is  very  complicated  and  beyond  the 
scope  of  this  report.  Here,  we  obtain  a  lower  bound  on  the  probability  distribution  of  suc¬ 
cessful  data  packets  by  pretending  that  all  packets  belong  to  different  data  nodes.  Clearly, 
since  any  intelligent  protocol  would  do  better  than  this.  (4.10)  serves  as  a  lower  bound  to 
the  true  probability  and  therefore  all  expressions  about  data  throughput  and  delay  based 
on  (4.10)  are  pessimistic.  For  moderately  small  values  of  the  data  node  activity  parameter 
Pd,  the  probability  of  having  more  than  one  packet  in  the  backlog  N?s  belonging  to  the 
same  data  node  is  relatively  small,  so  that  the  lower  bound  of  (4.10)  is  relatively  close 
to  the  actual  value.  For  larger  values  of  Pd  this  is  not  the  case  and  further  research  is 
required  to  determine  the  tightness  of  the  lower  bound  of  (4.10)  and  the  resulting  loss  in 
data  throughput  and  delay  performance. 

Equations  (4.4 )-( 4. 10 )  provide  us  with  the  transition  probability  matrix  of  the  two- 
dimensional  Markov  process  D. ' .  We  may  evaluate  the  stationary  distribution  of  the  process 
to  obtain  the  throughput -delay  performance  of  the  data  traffic  over  the  satellite  subnet. 

The  expected  data  throughput  over  the  S-subnet  is  then  given  by 

A*  Md  L-l-S'” 

i’d=  Y  Y  Y  5,‘e[s,'|A’.'",Ar]e[A’,'".A7']  (4.11) 

1VU*=0  V4._o  S'=  0 

»  I  * 


17 


in  packets  per  frame.  In  (4.11)  the  joint  distribution  P  is  obtained  from 


P  [N?s,N»a]  =  P  ( N?s  |  TV*'4)  P{N?') 


(4.12) 


where 


P  {Nfa\N?s)  =  ]T  P  {Nfa  |  C,_2iA74)  P  (C._2  |  A74) 


-  5]  p(a74  |c,_2)p(C,_2  IA74), 

C..7=Nf 


(4.13) 


and 


P{N?a  =J t  |  C,_2  =  m)  =  b(m,  k,~),) 


P(C,_2  I  A74)  = 


P(C,_2.A74) 


P(A74) 


ES-0  P(  c, ,  .V.“ ,  c,.2,A’"2  ) 

P(A\"*> 

E“=oEnI-,=<,P(C,.AV  I  Ci-2.A’~J)./>(C..,.A-~J) 

P(N?a) 


(4.14) 


(4.15) 


is  obtained  from  the  transition  probability  matrix  and  the  stationary  distribution  of  P4  = 
(C, ,  A7* )•  Moreover,  the  conditional  probability  P  (S4  |  Nfa,Nl'a)  is  obtained  from 


A/d  4*  m 


P  [5;  =  n  |  A',*14  =  m.  A’,*'4  =  *■]  =  ]T  P  [54  =  n  |  Nfa  =  m.  .V,1'4  =  *.  }',4  =  /] 

/=max { ni.n } 

•  P  [r,4  =  /  I  \fa  =  m.  A74  =  k]  (4.1G ) 


where 


p  [5;  =  77 1  xfa  =  m,  A74  =  k.y;  =  /]  =  p  [5, 4  =  n  |  A74  =  it.  r,4  =  /] 


is  independent  of  the  value  of  A',  and  is  given  by  (4.9)-(4.10).  and 

P  [r,4  =  /  |  \fa  =  m.  A* 4  =  Jt]  =  P  (.4;  =  /  -  m)  =  b(Md.l  -  m,Pd) 


(4.17) 


IS 


Finally,  the  average  data  delay  for  a  packet  successfully  transmitted  over  the  S-subnek. 
(measured  in  frames)  is  approximated  as 


D3  %  1  + 


E  {number  of  transmitted  data  packets  per  frame) 

E{ number  of  successfully  transmitted  packets  per  frame) 

1 


=  1  +  3  '  I™  E{Yn 


Vd 

ni 


=  1  +  p  '  ,‘i“  E  E  E  Y'p  lYwtM<K‘]  p 

J  '  °°  N,V*  Nf>  Y‘ 


hi  oc  Mi  + 


2  V  ^  d  T 

= 1  +  -j  •  ^  E  E  E  lp I1’.’  =  <i\'j  =  m.A’r  =  ‘-]  ^  [ay*  =  m,.vr  =  t] 

'  °°*=0m=0  /=r~ 


(4.18) 

where  P  [F,J  =  /|jV,dj  =  m^N?3  =  fc]  is  obtained  from  (4.17)  and  P  [Artds  =  m.  N^3  =  Jt]  is 
obtained  from  (4.12).  The  term  1  in  the  right  hand  side  of  (4.18)  accounts  for  the  fixed 
delay  of  one  frame  (=  L  slots)  introduced  by  the  satellite  to  all  data  packets  transmitted 
through  it. 

Notice  that  the  expression  in  (4.18)  is  only  an  approximation  to  the  actual  data  delay 
over  the  S- subnet.  The  method  followed  for  this  derivation  parallels  the  approximate  delay 
analysis  provided  in  [14,  Section  4.2,  pp.  219-221]  for  slotted  multiple-access  systems  with 
retransmissions.  Little’s  law  is  not  exactly  used  in  (4.18),  but  rather  a  similar  intuitive 
argument  produces  the  first  equation  in  (4.18):  refer  to  [14]  for  a  detailed  description  of 
this  approximation.  Our  situation  is  different  from  that  of  [14]  in  that  we  are  interested 
in  the  overall  system  delay  (waiting  time  +  service  time)  whereas  [14]  deals  with  waiting 
time  (queuing  delay)  only;  of  course  the  protocol  of  [14]  is  also  different. 

4.2.  Data  Traffic  over  the  G-Subnet 

Some  of  the  data  nodes  that  were  involved  in  a  collision  even  over  the  S-subnet  choose 
to  retransmit  the  collided  packets  over  the  ground  subnet. 

Let 
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R dx9  =  Number  of  data  packets  in  backlog  over  G-subnet  at  the  beginning  of  frame  i 
N*9  =  Number  of  data  packets  from  S-subnet  that  join  backlog  of  data  over  G-subnet 
at  the  beginning  of  frame  i 

S9  =  Number  of  successfully  transmitted  data  packets  over  G-subnet  in  frame  z 
T,  =  Number  of  packets  that  actually  transmit  over  G-subnet  in  frame  i 
Figure  5  depicts  all  components  of  data  traffic  over  the  ground  subnet. 

It  is  easy  to  see  that 

J??*,  =  Rdt9  +  X?9  -  Sf  (4.19) 

The  protocol  observed  by  the  data  terminals  over  the  G-subnet  is  as  follows:  At  the  begin¬ 
ning  of  frame  z,  the  terminals  with  packets  to  transmit  in  the  G-subnet  have  information 
about  the  number  of  established  voice  calls  and  the  data  backlog  over  the  G-subnet.  The 
total  number  of  these  data  terminals  is  R dt9  4-  Xf9\  each  of  these  data  nodes  transmits  its 
packet  in  frame  i  with  probability  /,  given  by 


lO;  if  < 


(4.20) 


where 

L'(N;’9)  =  K9d  -  X'9  (4.21) 

denotes  the  unused  multiple-access  capability  of  the  CDMA  channel  of  the  G-subnet  in 
frame  i:  this  channel  capacity  is  available  for  use  by  the  T,  (out  of  the  R f9  -I-  X,9)  data 
users  that  actually  transmit. 

The  selection  of  /,  in  (4.20)  as  the  retransmission  probability  is  well  justified  by  the 
fact  that  it  makes  the  average  number  of  transmitting  users  equal  to  the  available  CDMA 
capacity,  that  is. 


E{T,}  =  f,(Rdx9  +  X?9 )  =  L'(X?9)  =  K9  -  X;6 


This  type  of  retransmission  probabilities  were  first  introduced  in  [10]  and  [11]  in  the  context 
of  ALOHA  protocols  with  retransmission  control.  We  have  also  used  them  in  [G]  for  the 
retransmission  control  of  multi-media  (voice/data)  CDMA  networks.  Actually,  several 
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variants  of  (4.20)  were  used  together  with  different  feedback  information  for  the  schemes 
termed  A,  B,  and  C  whose  performance  was  analyzed  in  [6]. 

Let  D9  =  (Kf9,Nf9,N*9,Ci-i).  Then 
P\H9t+i  ~  (rni,m2,rn3,mi)\D9  =  (ki,k2,k3tk4)] 

=  P[R%,  =  |  Nf*x  =  m2,  A7"9,  =  m3,C,  =  m4,.R?9  =  A N*9  =  A-2, N*9  =  fc3,C,_,  =  *4] 

•  I  Ng ,  =  m3,C.  =  m4,i?f9  =  A-, , A7,**9  =  &2,  A7”9  =  *3,C,-_,  =  k4) 

■  W+’i  =m3\C,  =  m4,Rd9  =  fc,,^9  =  k2,N?9  =  k3,C,-x  =  M 

■  P[Ci  =  m4  |  Rdx9  =  kuNd9  =  A:2 ,  Nx9  =  *3lC._,  =  *4]  (4.22) 

For  the  third  of  the  four  terms  in  the  sum  on  the  right  hand  side  of  (4.22),  since  voice 
evolves  independently  of  data,  we  have: 

*W+i  =  m3  I  C,  =  mA,Rd9  =  ki ,  Nx9  =  k2,  A',*9  =  h,Ct-i  =  M  =  P[N^X  =  m3  |  A’,1’9  =  k3], 

(4.23  a) 

which  is  obtained  from  the  joint  conditional  distribution  P{ N^1 ,  A7*^  IA7,1'*,  A7,1’9)  and  the 
stationary  distributions  P(Nxa,Nx9)  and  P(N x9 )  evaluated  in  Section  3  as 


p(n:’,in:’)  = 


P(N ’•) 

V  O/  \-rs  V*’9  \’t’S  \Tl’9  \ 

Nv'=0  2-jNv,=0  ■r'-'i+l'*'i+l>*Ni  1 


■  +  i 


F(.v;y) 


E!T”=i  £&.,  e(.v, 7.  IA?*.Ar)  •  mr-NT) 

I  -f  1  » 


Ztt.=0r<xr<K9) 


(4.23b) 


Moreover,  for  the  second  term  on  the  right  hand  side  of  (4.22)  we  have 
/’[A'J,  =  m2  |  A',7,  =  m3.C,  =  nu.R*9  =  k^.Xf9  =  k2 .  A'*’9  =  *3.C.-,  =  *4] 

=  ^iArtd+9,  =  m2  |  Ct_,  =  A-.,] 

=  b(k4,m2, 1  -  t-.+  i  ).  (4.24) 

In  deriving  (4.24)  we  used  the  fact  that  that  A’^.  the  number  of  data  packets  from  the 
satellite  subnet  joining  the  backlog  of  the  ground  subnet  during  frame  i  +  1,  depends  only 
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on  C,_  i,  the  number  of  colliding  data  packet  over  the  satellite  subnet  during  frame  t  —  1; 
refer  to  (4.3)  for  the  expression  determining  the  splitting  probability  7,+i- 
The  fourth  term  on  the  right-hand  side  of  (4.22)  is  given  by 

P[C ,  =  m4  |  Rdg  =  kuNdg  =  k2 , N"9  =  k3,C,-i  =  k4] 

hi  OO 

=  T  V  P[Cl  =  m4\Nr,Cl^Rd3  =  h,Ndg  =h,N”0  =  h,C,-i=k4] 
^"*  =  0  c,_2- o 

•  P[C,_2  |  N,s ,  Rdg  =  k,,Ndg  =  k2,N"9  =  k2,  C,—i  =  k4 } 

•  P[N r  I  Rdg  =  kuNda  =  k2,N"9  =  h,C,-i  =  k4] 

h'l  oo 

=  T  V  P{C,  =  m4  |  A7S. C, _2]  •  P\Ct-2  |  N dg  =  k2\  ■  P[Nf  |  N’9  =  h) 


(4.25a) 


where  the  first  term  on  the  right  hand  side  of  (4.25a)  can  be  obtained  from  (4.5),  the 
second  term  can  be  obtained  from  (4.24)  and  is  given  below  as 


P(C,_2  |  Nr,Rl9,Nd3,N”9,C,-i) 

^  -P(  C. - 2  ^7 J ,  Rd,9 .  N?9 .  N*9 .  C,  _  i ) 

P(  N** ,  Pf 9 ,  Ndg ,  N*9 ,  C,  _  i ) 

=  P( A7,J. Pf9,  A’,*9,  AVg.  C,-,  (C,-2 )  •  P(C,-2 ) 

P(Nl’\Rd9,Nd9.N,,gX,,-l) 

P( Ndg  1 C, _2 .  A7a •  Pf g •  A79 .  C, _ j )  •  P(A7a.Pf9.A7g,C,-ilC,-2)  •  P(C,-; 

P(A7s.Pf9.A'Irf9.A7?.C,_1) 


P(  A^9 1 C,  - 2 )  -  P(  .v,"\  pf 9 .  A',"g .  C,  - 1 1 C’,  -  2 )  ■  P(  C,  _ 2 ) 
P(  A’7 .  pf 9 .  A*,"9 .  C,_ ,  |  A';'9 )  •  P(  Ndg ) 

P(Nd9\C,-2 )  •  P(A7*.Pf9.  A79.C',-1 )  •  P(C,_2) 
P(A7*.Pf9.A79.CI_,  )  ■  P(A',d9) 

PfA^9  I  C,_2)-P(C,.2) 

P(  Ar,<<9 ) 
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=  P(C,_ 2  I  N*9) 

P(Nf9  1  C,_2)-P(C,-2) 
E“.J=0  P(^s|c,_2)-P(c,_2) 


fe(C,-2,A^l-7,)-P(C,-2) 

^  i 


(4.256) 


where  P(C,_2)  is  obtained  from  the  Markov  chain  ( Ct,N *'J)  of  Section  4.1,  and  the  last 
term  on  the  right-hand  side  of  (4.25a)  is  given  as 


P(N? 


N-9 


)  = 


p(n?\n; 


vg  , 


P(N?\N"9) 


p(n; 


~vg 


(4.25c) 


Finally,  for  the  first  term  in  the  sum  on  the  right  hand  side  of  (4.22).  we  have  two 
options  for  the  model  of  correct  reception  of  data  packets:  the  threshold  model  and  the 
graceful  degradation  model.  If  the  threshold  model  of  Section  2.1  is  used  for  the  reception 
of  data  packets,  according  to  which  data  packets  are  successful  (i.e.,  correctly  received)  if 
the  total  number  of  simultaneously  transmitted  packets  (data  and  voice)  over  the  G-subnet 
is  smaller  than  A'J,  we  have 


P[Rd,U  =  m,  |  K&  =  =  m3.C,  =  m4.Eft9  =  X?9  =  k2,X,"9  =  k3.C,-i  =  k4) 

-P\Rdi+i  =  =  ^'i  •  x*9  =  k2'N't9  =  k3] 

=  1  =  ™ilT-  =  l’Rd,9  =  h'N?9  =  h.x;9  =  k3)  •  p\t,  =  i\Rdt9  =  ki.xf9  =  k2.x;'9  =  k3\ 


i 


-!- 6( kx  -(- (),  /, )  if  mi  =  A  ]  +  A-2.  0  <  £/(  A’3  )  <  k]  +  A-2 

1  if  77!  J  =  A'l  A'2  .  L'  (  k3  )  <  0 

6(  k\  -(-  A’2 ,  k'i  k~2  —  m  i .  /, )  if  max {0,  k\  4-  k2  —  L1  ( k3  )}  <  m  ]  <  A-]  -+•  .I'2 

I  0  elsewhere 


where  L'{k3 )  =  A’J  -  A'3  as  defined  by  (4.21 ). 


(4.26a) 
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On  the  other  hand,  we  may  use  the  graceful  degradation  model  for  CDMA 
multiple-reception  of  data  packets  (see  [10]),  according  to  which  there  is  a  non-zero  proba¬ 
bility  P(n\k)  of  correct  reception  for  any  arbitrary  number  n  of  packets  and  k  is  the  total 
number  of  simultaneously  transmitted  packets  (data  and  voice)  over  the  G-subnet.  In  this 
case  we  have 


**[*?+ 1  =  |  AT&  =  m2.A^,  =  m3,C,  =  m4,Rd9  =  ku  N dg  =  k2.N't'9  =  *3,C\- ,  =  k4] 

=  P\Rlh  =  mAR*9  =  kuN?9  =  k2,N-'9  =  k3) 

=  Y^P[R‘!U=m1\Tt  =  LRdtg  =  kl,  Ndg  =  A' 2 ,  =  *3]  -  P[T,  =  l\Rdg  =  .  JV^  =  A-.,  Ar,r?  =  jt3] 


—  P  (A-j  +A'2  —  mi  successful  data  packets] /-fA’3  total  transmitted  packets)  6(  A'j  +k2.  l.f,  ] 


1=0 


(4.20  b) 


where 


P  (n  successful  data  packets  |  /  +  k2  total  transmitted  packets) 

^  j  (;)[1_P£(/+ A'3)]n[P£(/  +  A-3)]'-n  if  0  <  n  <  / 
'  0  elsewhere 


(4.27) 


and  the  probability  of  packet  error  for  a  frequency-hopped  CDMA  system  using  Reed- 
Solomon  (A.  A  )  coding  and  erasures-decoding  is  given  for  k  simultaneous  packet  trans¬ 
missions  by 

PE(k}=  Y,  (*)[Ph(k)Y[l-P,.(k)]s-'  (4.2S) 

i=  .V  —  h’+  1  '  ‘  ' 

where  the  probability  of  a  hit  from  other  users  is  given  by  [11] 


Pk(k) 


1  - 

1  - 


0-0 

0-0 


k- 1 

A-l 


for  hoj)  synchronous  FH  systems 
for  hop  asynchronous  FH  systems 


(4.29) 
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The  approximation  in  (4.27),  termed  the  independence  of  receiver  operation  assump¬ 
tion  (IROA),  has  been  shown  in  [10]  to  have  satisfactory  accuracy  for  most  situations  of 
practical  interest.  The  exact  expression  for  the  quantity  in  the  left  hand  side  of  (4.27)  has 
computation  complexity  that  is  typically  exponential  in  n  (e.g.,  2"  —  1  sums  are  involved  in 
the  exact  expression  for  P(n|/)  for  FH/SS  systems  with  RS  coding  and  erasures-decoding); 
refer  to  [10]  for  a  detailed  discussion  of  these  issues  and  numerical  examples  establishing 
the  accuracy  of  the  approximation  of  (4.27).  Finally,  recall  that  our  system  is  packet  (  = 
frame)  synchronous  in  the  ground  subnet  and  slot  synchronous  in  the  satellite  subnet; 
however,  it  is  not  required  that  the  FH/SS  system  employed  in  the  ground  subnet  is  syn¬ 
chronous  at  the  hop  (or  dwell  time)  level;  this  is  the  reason  for  providing  two  expressions 
for  the  probability  of  a  hit  in  (4.29). 

Finally,  notice  that  the  number  of  data  packets  in  backlog  at  the  beginning  of  frame 
i  over  the  G-subnet,  i?f9,  has  a  countably  infinite  space;  since  it  is  possible  (although 
with  very  small  probability)  that  several  of  the  data  packets  that  get  transmitted  with 
probability  /,  are  not  actually  transmitted  or  they  collide  repeatedly  and  thus  according 
to  our  packet  arrival  model  data  packets  may  accumulate  indefinitely  in  the  buffers  of 
the  data  nodes.  Therefore,  in  the  numerical  evaluation  we  must  truncate  the  transition 
probability  matrix  of  (4.22)  and  (4.26)  to  some  moderately  large  value  of  R^9 . 

The  expected  data  throughput  over  the  G-subnet  can  be  computed  as 


h  d 


L‘ { .\'c 


=  y  Y  Y  Y  s9 P( s>  i Rd*' 9 ’ N>9 • N'9 )P[ R>9 ■  Y9 ■  -vr ) 


(4.30) 


.\v,=0  R^o-0  Sf=  0 


in  packets/frame  and  for  the  threshold  model  of  CDMA  other-user  interference 


P(  S f  =  n\Rdtg  =  *, .  X?9  =  k2 .  A',*’9  =  k3  ) 

=  Y  P(S?  =  n\T<  =  l<  K i9  =  ti  •  Y9  =  Ks  =  h)-  P(T,  =  l\R*9  =  .  Y9  =  k2.  A’."9  -  h 
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£f=l‘*a)-H  +  *2,  lji)\  +b(ki  +k2,0ji)  if  n  =  0,  0<L'(k3)<k1+k2 

1  if  n  =  0,  L'(Jfc3)<  0 

b(kx  +  k2,n,f,)  if  0  <  n  <  min{L'(A:3),  h  +  M 

0  elsewhere 

(4.31a) 

where  L'(k3)  is  defined  as  in  (4.21).  If  the  graceful  degradation  model  for  CDMA 
other-user  interference  is  employed  we  should  use 

P(S?  =  n\Rdl9  =  k\ ,  Ndg  =  k2,  N"9  =  k3) 

=  £P(S?  =  n|Ti  =  /,P?9  =  A,,^9  =  A-2.iV;’s  =  k3) 

1 

•  P(Tt  =  /|Pf9  =  Jc,,Arfs  =  A2.:V;9  -  *3) 

*1  +*2 

=  P(n  successful  data  packets  |  /  +  A-3  total  transmitted  packets)  b(kt  -f  A*2,/, /,) 
1=0 

(4.31 b) 

where  P(n  successful  data  packets  |  l  k3  total  transmitted  packets)  is  given  by  (4.27)- 
(4.29).  Finally,  on  the  right  hand  side  of  (4.30)  P{Rdg  =  J.X?9  =  9  =  n)  denotes 

the  stationary  distribution  of  the  Markov  chain  D9  whose  one-step  transition  probabilities 
were  evaluated  above. 

Finally,  the  expected  data  delay  of  a  packet  successfully  transmitted  over  the  G-subnet 
(measured  in  frames)  is  given  by 

D9  =  D9*  +  D"  (4.32) 

where  D 99  and  D"  denote  the  components  of  the  delay  of  a  typical  packet  transmitted 
over  the  G-subnet,  which  are  due  to  protocol  operation  over  the  S-subnet  and  G-subnet. 
respectively,  and  are  given  by 
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and 


D"  ml + 


E  {number  of  backlogged  data  packets  per  frame  in  G  —  subnet) 
E{ number  of  successful  data  packets  per  frame  in  G  —  subnet) 
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In  (4.33).  P  [-V,'’’’  =  m.Nf"  =  k]  and  P[C,  =  n\.\}'!'  =  ni.Y *  =  /]  are  obtained  from 
(4.12)  and  (4.9)-(4.10).  respectively.  Similarly,  in  (4.34).  P(Rdt9 ,Rdg )  can  be  evaluated 
from  P(Rdt9.  N?9,  N't9  ),  the  stationary  distribution  of  the  Markov  chain  JD9  =  (i?f9.  X?9 .  N't  9  ). 
Moreover,  the  term  1  on  the  right  hand  side  of  (4.33)  accounts  for  the  fixed  satellite  delay 
of  one  frame  that  any  data  packet  transmitted  over  the  satellite  is  subject  to.  The  second 
term  in  (4.33)  represents  the  component  of  the  delay  due  to  collisions  and  retransmissions 
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over  the  S-subnet  before  a  data  packet  gets  retransmitted  over  the  G-subnet.  The  term  1 
in  the  right  hand  side  of  (4.34)  accounts  for  the  fact  that  data  packets  from  the  S-subnet 
retransmitted  through  the  G-subnet  are  delayed  by  an  additional  one  frame  period  [because 
of  the  delay  in  the  satellite  channel  feedback  providing  information  about  packet  collisions]. 
Finally,  the  second  term  in  the  right  hand  side  of  (4.34)  represents  the  component  of  the 
delay  due  to  collisions  and  retransmissions  over  the  G-subnet. 

Finally,  notice  that  the  expressions  in  (4.33)  and  (4.34)  are  only  approximations  to 
the  actual  delay  terms.  The  method  followed  for  their  derivations  is  similar  to  that  used 
for  equation  (4.18)  (data  delay  over  the  S-subnet)  and  parallels  the  approximate  delay 
analysis  provided  in  [14,  Section  4.2,  pp.  219-221]  for  slotted  multiple-access  systems  with 
retransmissions. 

5.  NUMERICAL  RESULTS 

In  all  figures  with  numerical  results  presented  in  this  section,  the  parameters  of  voice 
and  data  traffic,  as  well  as  the  capacities  of  the  ground  and  satellite  channels,  are  depicted 
on  each  figure. 

The  offered  voice  traffic  Gv  is  defined  as 

A'f+A'i 

Gv  =  lim  V  [(Mv  -  n)Pv  +  n]P  [N*9  +  N*s  =  n] 

l  —  OO 

n=0 

Ki+h’l 

=  MVPV  +  (1  -  Pv)  y  n  lim  P  \N;'9  +  Nr  =  n]  (5.1) 

'  I— *OC 

n=:0 

wher  lim.-.oc  P{ N]  9  -I-  A7,13  =  n]  is  given  by  (3.12). 

The  basic  set  of  system  and  traffic  parameters  that  is  used  in  all  figures  presented  below 
is:  Mv  =  20  and  Mj  =  15  for  the  total  population  of  voice  and  data  users,  respectively: 
K9  =  7  and  I\9d  —  5  for  the  multiple-access  capability  of  the  CDMA  ground  network  [refer 
to  (2.1 )-( 2.2 )]  for  voice  and  data  users,  respectively;  p  =  0.05  for  the  parameter  of  the 
geometric  distribution  governing  the  length  of  the  voice  calls  [refer  to  (2.3)]:  and  L  =  13 
for  the  frame  length,  which  is  equal  to  the  constant  delay  introduced  dv  the  satellite. 

Voice  throughput,  defined  as  the  expected  total  number  of  voice  calls  in  progress  per 
frame  (summed  over  the  G-  and  S-subnets),  is  shown  in  Figure  6  as  a  function  of  the  offered 
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voice  traffic  Gv  defined  by  (5.1).  All  voice  calls  jure  accepted  when  the  voice  traffic  load 
is  light.  For  a  heavy  voice  traffic  load,  the  throughput  reaches  saturation  as  the  number 
of  new  voice  calls  increases,  since  most  of  the  new  calls  are  blocked.  The  traffic  level  at 
which  saturation  is  reached  depends  on  A'*,  the  number  of  slots  per  frame  in  the  satellite 
channel  that  are  available  for  voice  traffic.  Clearly,  the  availability  of  time  slots  in  the 
satellite  channel  increases  the  voice  throughput  of  the  overall  system. 

Figure  7  illustrates  the  blocking  probability  of  voice  calls  in  the  satellite  and  ground 
networks,  as  a  function  of  the  offered  voice  traffic  Gv.  As  we  allocate  a  larger  part  of 
the  satellite  channel  capacity  to  voice  traffic  (by  increasing  A'*),  blocking  performance 
improves.  However,  an  increase  in  A\*  results  in  a  decrease  in  the  number  of  slots  available 
for  data  traffic  in  the  satellite  subnet,  and  thus  in  the  degradation  of  the  data  performance 
measures. 

Average  data  throughput  over  the  satellite  subnet  (S-subnet)  is  shown  in  Figure  8 
as  a  function  of  the  duty  factor  Pd  of  the  data  nodes.  In  Figure  8.  when  voice  traffic  is 
light,  as  is  for  the  case  ( Gv  =  3.6,  Pv  =  0.01,  A'*  =  5),  there  are  more  channels  available 
for  data  traffic,  fewer  collisions  occur  and  throughput  over  the  S-subnet  initially  increases 
and  then  starts  decreasing  slowly  as  increases  beyond  the  value  0.35.  due  to  increasing 
data  packet  collisions.  On  the  other  hand,  when  voice  traffic  is  heavier,  as  is  for  the  case 
(G„  =  11.7,  Pv  =  0.1,  A'*  =  5),  data  throughput,  after  reaching  its  (significantly  lower) 
maximum  value,  decreases  due  to  an  increase  in  the  number  of  colliding  packets.  Finally, 
for  the  case  of  substantially  heavier  voice  traffic  ( Gv  =  18,  Pv  —  0.5.  A'*  =  10),  there  are  so 
few  channels  available  for  data  traffic  that  the  data  throughput  decreases  slowly  over  the 
entire  range  of  values  of  Pd-  Changing  the  routing  parameters  a  and  3  [refer  to  equation 
(4.3)]  does  not  significantly  affect  the  data  throughput  for  either  light  or  heavy  voice  traffic. 
It  appears  that  for  light  voice  traffic  (G,.  =  3.6)  there  is  some  loss  in  the  data  throughput 
if  values  of  a  and  3  that  divert  more  data  traffic  to  the  G-subnet  (e.g.,  o  =  1.0  and  3  =  3.0 
versus  a  =  0.6  and  3  =  8.0)  are  used:  this  observation  is  valid  for  the  entire  range  of  values 
of  Pd-  As  the  offered  voice  traffic  load  increases  to  G,.  =  11.7  the  pair  (a,  3)  =  (0.6. 8.0) 
results  in  slightly  larger  data  throughput  than  the  pair  (a,  3)  =  (1.0, 3.0)  for  values  of 
Pd  smaller  than  0.35,  whereas  the  opposite  is  true  for  larger  values  of  Pd-  Finally,  for 
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heavy  voice  traffic  ( Gv  =  18.0),  the  pair  {a,0)  =  (1.0, 3.0)  results  in  slightly  larger  data 
throughput  than  the  pair  (a,  ft)  =  (0.6, 8.0)  for  all  values  of  Pd\  that  is,  it  is  advantageous 
to  divert  more  data  traffic  to  the  G-subnet. 

The  parameters  a  and  0  of  (4.3)  determine  the  value  of  the  splitting  factor  1  —  7, +2, 
that  is,  the  fraction  of  the  data  traffic  of  the  S-subnet  involved  in  packet  collisions  (i.e.,  C, ) 
which  is  diverted  to  the  G-subnet.  From  these  two  parameters,  o  determines  the  threshold 
(which  is  equal  to  a  fraction  of  the  overall  number  of  data  users)  which  should  be  surpassed 
by  C,  for  data  traffic  to  be  diverted  to  the  G-subnet,  whereas  0  determines  the  actual  value 
of  the  splitting  factor  as  /?/C,.  The  effect  of  the  splitting  factor  becomes  more  pronounced 
as  a  and  0  increase  because  in  this  case  both  the  threshold  and  the  value  of  the  splitting 
factor  increase.  In  most  situations,  larger  values  of  a  cause  larger  increases  in  1  —  7, +2 
than  larger  values  of  3. 

Figure  9  displays  the  average  data  delay  over  the  S-subnet  as  a  function  of  the  duty 
factor  of  the  data  nodes  for  the  same  system  and  traffic  parameters  as  in  Figure  8.  Similar 
trends  as  in  Figure  8  are  observed.  For  light  voice  traffic  the  data  delay  is  small  and  is 
almost  insensitive  to  the  changes  in  the  values  of  o  and  0.  However,  for  heavier  voice 
traffic  the  data  delay  increases  and  becomes  more  dependent  on  the  values  of  a  and  0\ 
actually,  as  more  data  traffic  is  diverted  through  the  G-subnet  [(a,/?)  =  (1.0. 3.0)  versus 
(a,0)  =  (0. 0,8.0)]  the  data  delay  over  the  S-subnet  improves  significantly. 

In  Figures  10  and  11,  the  data  throughput  and  data  delay  over  the  ground  subnet 
(G-subnet)  are  illustrated  as  functions  of  the  duty  factor  of  the  data  nodes  Pd  for  the 
thresho'  '  model  described  in  Section  2.1.  The  parameters  of  the  threshold  model  for  the 
ground  s,  net  are  A'-?  =  7  and  I\9d  =  5.  These  numbers  correspond  to  the  multiple-access 
capability  of  a  specific  FH/SS  system  (with  parameters  to  be  provided  below)  for  desirable 
packet  error  probabilities  of  approximately  =  10-3  (for  voice)  and  P g  =  10--’  (for 
data)  in  (2.1)  and  (2.2),  respectively.  The  aforementioned  FH/SS  system  employs  32- 
ary  (A/  =  32)  FSI\  modulation  with  noncoherent  demodulation,  40  hopping  frequencies 
(i.e.,  q  =  40  in  (4.29)].  and  an  RS  (32.16)  code  [i.e.,  N  =  32,  A  =  16  in  (4.28)]  with 
erasures-decoding;  perfect  side  information  is  assumed  available  from  the  channel,  that 
is,  all  erasures  are  only  due  to  multiple-access  interference  and  can  be  detected  without 
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errors.  The  other  system  and  traffic  parameters  are  the  same  as  the  ones  used  in  Figure  8. 

In  Figure  10  in  which  the  data  throughput  over  the  G-subnet  is  plotted  versus  Pj, 
when  voice  traffic  is  light,  the  number  of  channels  available  for  data  increases  so  that 
throughput  increases  as  Pd  increases  to  0.3  and  then  decreases  slightly.  On  the  other 
hand,  when  voice  traffic  is  heavier,  data  traffic  can  use  only  the  channels  that  axe  not 
used  for  voice  calls  in  the  G-subnet.  Data  throughput  thus  decreases  faster,  because  the 
number  of  colliding  data  packets  increases  whenever  voice  traffic  (Gv )  and  data  traffic 
(Pd)  increase.  Comparison  of  the  throughput  for  two  different  pairs  of  the  parameters 
(a,0)  shows  that  the  pair  (a,/3)  =  (0.6, 8.0)  results  in  better  performance  than  the  pair 
(1.0, 3.0)  suggesting  that  an  overall  smaller  splitting  factor  (1  —7,)  is  preferable  because 
fewer  data  packet  collisions  occur  in  the  G-subnet  (due  to  fewer  packets  being  diverted 
from  the  S-subnet)  and  more  data  traffic  goes  through  successfully. 

In  Figure  11,  the  average  data  delay  over  the  G-subnet  (which  includes  the  delay 
experienced  while  unsuccessfully  attempting  transmission  over  the  satellite  subnet)  is  dis¬ 
played  as  a  function  of  the  duty  factor  of  the  data  nodes.  The  same  threshold  model  as 
for  Figure  10  is  used.  In  general,  data  delay  increases  as  voice  traffic  becomes  heavier  over 
the  G-subnet.  The  splitting  factor  has  a  significant  effect  only  when  the  voice  traffic  is 
heavier,  in  which  case  larger  splitting  factors  cause  an  increase  in  data  delay. 

In  Figures  12  and  13,  the  average  data  throughput  and  data  delay  over  the  G-subnet 
are  shown  as  functions  of  the  duty  factor  of  the  data  nodes  for  the  graceful  degradation 
model.  This  model  is  described  in  detail  in  Section  4.2  [see  equations  (4.27 )-( 4.29 )] .  The 
same  FH/SS  system  model  parameters  used  for  the  threshold  model  (refer  to  the  discussion 
about  Figures  10  and  11  above)  are  used  here  as  well.  However,  to  reduce  the  computational 
burden  in  the  graceful  degradation  model  it  has  been  necessary  to  truncate  P(n  |  /  +  k 3 )  of 
(4.27).  In  particular,  in  the  process  of  generating  numerical  results  we  evaluate  P(n  j  / -f  k'3  ) 
from  ( 4.27 )-( 4.29 )  only  for  the  range  0  <  /  <  10  and  set  P{n  1  /  +  k'3 )  =  0  for  /  >  10  for  all 
0  <  fa  <  A'-?.  For  k'3  =  I\*  -  7.  this  corresponds  to  setting  the  packet  error  probability 
[of  (4.28)]  Pf(l  -+  I13)  =  1  for  all  values  of  PeU  +  k3)  larger  than  0.7.  For  k'3  =  0.  we  set 
PE{1  +  k3)  =  1  for  all  values  of  Pe(1  +  fa)  larger  than  0.0S.  For  each  intermediate  value 
of  k3  (0  <  fa  <  I\$)  the  truncation  (that  is  forcing  PeG  +  ^’3)  =  1)  takes  place  at  some 
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different  value  of  Pe(/  +  £3)  between  0.08  and  0.7.  Since  large  values  (close  to  1)  of  the 
packet  error  probabilities  Pe(/  +  t3)  result  in  extremely  small  probabilities  P(n  |  /  +  £3), 
truncating  the  graceful  degradation  model  in  this  way  does  not  significantly  reduce  its 
accuracy  for  high  voice  traffic,  because  in  this  case  k 3  takes  larger  values  (closer  to  A'*?) 
with  high  probability.  However,  for  low  voice  traffic,  in  which  case  A-3  takes  smaller  values 
with  high  probability,  the  accuracy  of  the  graceful  degradation  model  decreases  and  the 
results  of  any  comparison  with  the  threshold  model  are  expected  to  be  pessimistic. 

Since  the  graceful  degradation  model  described  above  is  computationally  very  de¬ 
manding,  in  Figures  12  and  13  we  present  results  for  only  two  values  of  the  average  offered 
voice  traffic  per  frame  Gv  (light:  Gv  =  3.6,  Pv  =  0.01,  A  *  =  5  and  heavy:  Gv  =  11.7, 
Pt.  =  0.1,  A'®  =  5);  the  values  of  a  and  0  are  1.0  and  3.0.  respectively.  It  is  observed  that 
the  behavior  of  the  data  throughput  and  delay  is  basically  similar  to  that  observed  in  Fig¬ 
ures  10  and  11  for  the  threshold  model  for  the  case  of  light  voice  traffic.  However,  for  the 
case  of  heavy  voice  traffic,  the  graceful  degradation  model  shows  more  drastic  decrease  in 
the  throughput  and  increase  in  the  delay  as  Pd  increases  beyond  the  value  .4.  Comparison 
of  the  two  curves  in  Figure  12  establishes  that  for  small  values  of  Pd  data  throughput  is 
greater  for  heavy  voice  traffic  than  for  light  voice  traffic.  This  is  justified  by  the  fact  that 
more  data  packets  use  the  G-subnet  in  the  former  case  than  in  the  latter  case  due  to  the 
increased  number  of  data  packets  diverted  from  the  S-subnet  to  the  G-subnet. 

In  Figures  14  and  15.  the  average  data  throughput  and  data  delay  over  the  G-subnet 
are  plotted  versus  the  duty  factor  of  the  data  nodes  for  both  the  threshold  and  the  graceful 
degradation  models.  These  comparisons  are  carried  out  for  low  (3.G)  and  high  (11.7) 
average  values  of  the  offered  voice  traffic  per  frame  Gv.  The  results  of  the  comparison  of 
data  throughput  in  Figure  14  are  different  for  the  low  and  high  G ,  values.  For  low  G,. . 
using  the  threshold  model  results  in  better  performance  than  using  the  graceful  degradation 
model,  for  all  values  of  Prf.  the  duty  factor  of  the  data  nodes.  This  can  be  explained  by  the 
fact  that  the  truncation  in  the  graceful  degradation  model  discussed  earlier  is  more  harmful 
to  it  for  low  voice  traffic.  If  such  truncation  were  not  necessary,  the  graceful  degradation 
model  would  perform  better  than  the  threshold  model  for  all  cases  of  interest.  By  contrast, 
for  high  G,..  using  the  graceful  degradation  model  results  in  better  performance  than  using 
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the  threshold  model  for  Pd  larger  than  0.25,  whereas  the  use  of  two  models  gives  very  similar 
results  for  Pd  smaller  than  0.25.  Moreover,  in  Figure  15,  the  data  delay  of  the  threshold 
model  appears  to  be  significantly  larger  than  the  delay  of  the  graceful  degradation  model 
for  both  low  and  high  Gv;  this  ordering  of  data  delays  is  valid  for  all  values  of  the  data 
duty  factor  Pd  except  for  values  smaller  than  .15,  in  which  case  the  delay  under  the  two 
models  is  nearly  the  same  for  light  voice  traffic. 

In  Figures  16  and  17,  the  average  data  delay  over  the  S-subnet  ( D 3 ),  the  average  delay 
over  the  G-subnet  (D9),  and  the  the  total  average  delay  of  a  typical  data  packet  in  the 
system  (transmitted  over  either  of  the  two  subnets)  D ,  defined  as 


D  = 


rid 


r)d  +  rd 


■  Da  + 


V9d 


rid  +  Vi 


■D9. 


(5.2) 


are  shown  as  functions  of  the  duty  factor  of  the  data  nodes  Pd  for  the  cases  of  light  voice 
traffic  (Gv  =  3.6,  Pv  —  0.01)  and  heavy  voice  traffic  (G„  =  11.7.  P,.  =  0.1).  respectively. 
The  various  delays  are  displayed  for  both  the  threshold  and  graceful  degradation  models. 
Clearly,  in  regard  to  comparisons,  the  total  average  delay  D  obtained  by  using  the  two  dif¬ 
ferent  models  (threshold  and  graceful  degradation)  of  correct  data  packet  reception  follows 
the  same  ordering  as  the  data  delays  over  the  G-subnet  discussed  in  Figure  15;  however, 
the  difference  in  the  values  of  the  total  delay  under  the  two  models  is  less  pronounced 
compared  to  that  of  the  delays  over  the  G-subnet. 


6.  CONCLUSIONS 

In  this  report,  we  have  studied  the  performance  of  certain  multiple-access  strategies  for 
voice/data  integration  in  heterogeneous  mixed-media  packet  radio  networks.  Our  models, 
protocols,  and  performance  analysis  are  applicable  to  networks  with  terrestrial  and  satellite 
links.  These  multiple-access  schemes  take  advantage  of  different  characteristics  of  the 
ground  and  satellite  radio  channels  to  better  serve  the  requirements  of  voice  and  data 
traffic.  In  the  ground  subnetwork.  CDMA  with  movable  boundary  in  the  code  domain 
is  used  to  serve  the  voice  traffic  and  a  fraction  of  the  retransmitted  data  traffic,  while 
framed  ALOHA  with  movable  boundary  is  used  on  the  satellite  subnetwork  for  data  and 
reservations  in  each  frame  are  used  for  voice.  In  our  scheme,  data  nodes  attempt  to 
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transmit  over  the  satellite  network  first  and  the  amount  of  data  traffic  transferred  to  the 
ground  network  depends  on  the  number  of  colliding  data  packets  over  the  satellite  network. 

An  accurate  Markovian  model  of  the  evolution  of  the  voice  and  data  traffic  over  the 
ground  and  satellite  subnetworks  has  been  developed  involving  two-  and  four-dimensional 
Markov  chains.  Moreover,  two  models  for  the  successful  reception  of  data  packets  are  used: 
a  threshold  model  and  a  graceful  degradation  model.  On  the  basis  of  the  above  models, 
the  steady-state  performance  of  the  two  subnetworks  is  evaluated  in  terms  of  the  blocking 
probability  of  voice  calls,  the  average  voice  throughput  (number  of  calls  in  progress  per 
frame),  the  average  data  throughput,  and  the  average  data  delay. 

The  performance  evaluation  of  the  multiple-access  schemes  considered  in  this  report 
has  established  the  interplay  between  different  traffic  types  at  different  offered  traffic  levels. 
It  has  also  established  the  sensitivity  of  data  throughput  and  delay  to  the  proper  selection 
of  the  parameters  determining  the  amount  of  retransmitted  data  traffic  routed  through  the 
ground  network;  this  is  more  pronounced  when  the  offered  voice  load  is  high.  Regarding 
the  threshold  and  graced.  '  .  gradation  models  considered  for  the  correct  reception  of  data 
packets,  we  have  show’,  chat,  for  low  offered  voice  traffic  loads,  using  the  threshold  model 
results  in  data  throughput  and  delay  performance  which  is  better  than  obtained  when  the 
graceful  degradation  model  is  used,  while  the  opposite  is  true  for  higher  voice  traffic  loads. 
If  such  truncation  were  not  nectary,  the  graceful  degradation  model  would  perform  better 
than  the  threshold  model  for  all  cases  of  interest. 

Overall,  our  modeling  and  performance  analysis  study  of  mixed  media  packet  radio 
networks  with  voice  and  data  traffic  can  provide  the  tools  for  understanding  these  networks, 
as  well  as  guidelines  for  the  proper  design  of  suitable  multiple-access  schemes. 
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G  -  Subnet 


D:  Active  Data  Nodes 
V:  Active  Voice  Nodes 
I:  Idle  Nodes 

_  Links  with  Direct  Transmissions  or  Retransmissions 

'  - - Links  with  Retransmissions  only 


Figure  1.  Mixed-Media  (Satellite/Ground)  Network 
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Figure  2.  Satellite  Subnetwork  Using  Framed  ALOHA 
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Figure  3.  Ground  Subnetwork  Using  CDMA 
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Figure  4.  Voice  Traffic  over  Ground  and  Satellite  Subnets 


39 


Average  Number  of 


Sf 


Figure  5.  Data  Traffic  over  Satellite  and  Ground  Subnets 
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MEAN  OFFERED  VOICE  TRAFFIC/FRAME  (Gv) 
Figure  6  Voice  throughput  over  S-  and  G-  nets 
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MEAN  OFFERED  VOICE  TRAFFIC/FRAME  (Gv) 
Figure  7  Voice  blocking  probability  over  S-  and  G-  nets 
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Figure  10  Data  throughput  over  G-net  (Threshold  Model) 
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DATA  NODE’S  DUTY  FACTOR  Pd 
Figure  14  Data  throughput  over  G-net  for  the  threshold 
and  graceful  degradation  models 
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